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Abstract 

We investigate the production of the weak bosons W~^, W^ and Z at the LHC as a function of their 

rapidity, reconstructed experimentally from their leptonic decay. We show that the measurements 

provide a powerful tool for constraining the parton distribution functions, as was already the case 

^5 ' for the lower energy pp collider data. We study the charge asymmetry determining the u and d 



Ph. 
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distribution functions using the reconstructed W~^, W rapidities. We also show how the ratio of 



^ . the W and Z boson rapidity distributions directly probes the charm quark distribution function 

CN ' of the proton. 
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I. INTRODUCTION 



The weak bosons W~^, W~ and Z are produced in a hadron collider predominantly 
through the interaction of a quark and antiquark from each hadron, the so-called Drell-Yan 
mechanism [l|. This process has been studied extensively as a test the Standard Model 
(SM) I2H17I. 

The cross sections for W^, Z production have been computed in quantum chromody- 
namics (QCD) up to the next-to-next-to leading order (NNLO) 18|, [lOj and provide a test of 
higher order QCD corrections. The measured cross sections are related to the partonic QCD 
cross sections by the parton distribution functions (PDFs). The partonic cross sections and 
the PDFs factorize making a direct study of the PDFs possible. 

Particularly, the PDFs can be examined using the rapidity dependence of the W and 
Z production cross sections. The method is familiar from the analysis of lower energy 
pp collider data 20| where it has been extensively used to investigate the u{x) and d{x) 
distributions using the ratio of the W^ and W~ rapidity distributions, referred to as the 
charge asymmetry 20l-|25l| . Additionally a prescription was proposed to determine the charm 
quark structure function from both W and Z rapidity distributions measured at the Tevatron 



26| . where the charm quark contribution to W, Z production is smal 



method, the strange quark distribution has been studied at the LHC 27 



Using a similar 



28 
201 



In this paper, we investigate the PDFs by applying the same methods [20|, |26| to W and 
Z rapidity distributions at the LHC center of mass energy ^/s = 7 TeV, to be increased to 
14 TeV in the future. These are significantly higher energies than the Tevatron energy of 
^/s = 1.8 TeV for Run I and 1.96 for Run II. The higher energy not only leads to increased 
production rates, it provides access to the lower fractional momenta of the partons where 
sea quark interactions dominate; these were difficult to probe at Tevatron energies. We 
will return to this later on. In addition to the increased energy, the LHC collides protons, 
a symmetric process, resulting in a different pattern of rapidity distributions and charge 
asymmetry compared to pp collisions. Specifically, using the charge asymmetry and the 
ratio of the W and Z boson rapidity distributions, we will show how the LHC data, unlike 
the lower energy Tevatron data, directly probe u{x)/d{x) and the charm quark distribution 
function, respectively. 



This paper is organized as follows. In the following section, we introduce the formalism. 



In Section III, we investigate the u{x), d{x) PDFs using the charge asymmetry of the exper- 
imentally reconstructed W~^ and W~ rapidity distributions. In Section IV, we demonstrate 
the sensitivity of the charm quark PDF to the ratio of the rapidity dependence of the W 
and Z cross section. Finally, in Section V, we summarize our results. 

II. RAPIDITY DISTRIBUTIONS OF THE DIFFERENTIAL CROSS SECTIONS 



The general expression for the cross sections for W, Z production is given by 

aAB^W/Z = ^ dXa dXbfa/A{Xa,Q'^)fb/B{Xb,Q'^)crab^W/Z , (l) 

a,b 

where the partonic scattering cross section is CTab^w/z- In Eq. [1], Xa{b) are the momentum 
fractions carried by the partons a{b) in the colliding hadrons A{B), and fa/A{xa,Q'^) and 
fb/B{xb, Q"^) are the PDFs of a and 6, respectively. The scale Q^ is set to be the mass squared 
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w/z- 



of the produced boson, Q^ = M? 
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FIG. 1: The total cross sections for the production of weak bosons, W^^' and Z in proton-proton 
collision. The separate contributions from the valence-sea (dashed) and the sea-sea (dot-dashed) 
interactions are also shown as function of energy. 

In Fi g. [H we have calculated the cross sections for inclusive W, Z production using the 
VRAP 19(1 program. The cross sections are evaluated to NNLO using the MSTW2008- 
NNLO set of PDFs 29|]. The total cross section for W^ production is different from that 
for W~ production while they are identical in pp collision. The separate contributions 
from the valence-sea and sea-sea quark interactions are also included. Since there is no 
contribution from valence-valence quark interactions for pp collisions, the role of the sea 
quarks is enhanced and, as energy increases, the sea quark contributions become increasingly 
important. They are responsible for the increased production of W~ and Z relative to W^ . 



In addition to the total cross section, we are here particularly interested in the differential 
cross sections as a function of rapidity that sample the momentum fractions of the partons 
and are therefore optimal for constraining the PDFs. Including the heavy quark flavors, the 
differential cross section for W and Z boson production in pp collisions are given by 

da .rr+v\ T^ 2-110 F 

-^{PP -> W+X) = Kw—y=-XiX2 X 

{\Vudf [u{xi)d{x2) + d{xi)u{x2)] + \Vusf [u{xi)s{x2) + s{xi)u{x2)] 
+ \Vcsf [c{xi)s{x2) + s{xi)c{x2)] + \Vcdf [c{xi)d{x2) + d{xi)c{x2)] 
+ \Vubf [u{xi)b{x2) + b{xi)u{x2)] + \Vcbf [c{xi)b{x2) + b{xi)c{x2)]} , 

^{pp ^ W^X) = Kw^-^x,X2 X (2) 

{\Vud\^ [u{xi)d{x2) + d{xi)u{x2)\ + \Vus\^ [u{xi)s{x2) + s{xi)u{x2)\ 

\Vcsf' [c{^i)s{x2) + s{xi)c{x2)] + \Vcd\^ [c{xi)d{x2) + d{xi)c{x2)] 

+ \Vub\^ [u{xi)b{x2) + b{xi)u{x2)] + \Vcb\^ [c{xi)b{x2) + b{xi)c{x2)]} , 

da „-,- ^. 27iGf 

— {pp -, ZX) = Kz^^x,X2 X 

{gl[u{xi)u{x2) + u{xi)u{x2) + c{xi)c{x2) + +c{xi)c{x2)] 
+g1[d{xi)d{x2) + d{xi)d{x2) + s{xi)s{x2) + s{xi)s{x2) (3) 

+b{xi)b{x2) + b{xi)b{x2)]} , 

where 

gl = (l-8sin2^vi//3 + 32sin^^H//9)/2 , 

gl = (l-4sin2^vK/3 + 8sin^^H//9)/2 . (4) 

Here Gp is the Fermi coupling constant, and \Vij\ are the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements. The factor K^z is the ratio of the next-to-leadin g or der (NLO) to 



the leading order (LO) cross sections given by Kw,z — 1 + 87ras(M^^)/9 [30|. The parton 
momentum fractions, Xi and X2 are related to rapidity Y by 

Mw,Z Y Mw,Z _y f-s 

xi = — -=- e , X2 = — 1=- e , (5) 

where Mw,z are the masses of the W and Z bosons. Thus, Y covers the range —ln{^/s/M^^'l) 
to ln(A/s/Mw,z). For numerical evaluation, we use Mw = 80.399 GeV, Mz = 91.188 GeV, 



Gf = 1.166 xlO-5GeV"^ and sin'^ Ow = 0.23 |31| 



For the evaluation of the rapidity distributions we use the NLO cross sections supple- 
mented by the /^-factor instead of evaluating the NNLO directly. As mentioned in the 
introduction, the PDFs and the partonic cross sections factorize in their dependence on the 
measured hadronic cross sections; our focus is on the constraints on the PDFs from the ratio 
of the differential cross sections. The charge asymmetry and the ratio of d(jw±/daz are 
insensitive to the QCD corrections, therefore giving essentially the same results for different 
orders. 
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FIG. 2: The rapidity distribution of the differential cross sections for W and Z boson production. 
The left and right figures are for energies of y/s = 7 TeV and y/s = 14 TeV, respectively. The 



results are presented for CTEQ6.6(solid) 



NNLO (dotted) 



3^, MSTW2008-NLO (dashed) [29|, and MSTW2008- 



291] PDFs. 



The distribution of the differential cross sections in rapidity, evaluated following Eqs. 
EHSl are shown in Fig. [2] for both the CTEQ6.6 [32J and MSTW2008 [29] PDFs, and for 
two different center-of-mass energy a/s = 7 and 14 TeV. The differential cross sections 
are symmetric with respect to the sign of the rapidity with daw±{X) = daw±{—y), while 
daw-{y) = daw+{—y) in Pp collisions. The accessible rapidity range is \Y\ < 4.3 at i/i 
= 7 TeV, and \Y\ < 5 at y/s = 14 TeV allowing for the exploration of parton momentum 
fractions x as low as 1.7 xlO^^ and 4.2 xlO^^, respectively. The range of the rapidity and 
momentum fractions is much broader than for the Tevatron whose rapidity range from —3 
to 3 probed minimum values of xi and X2 of 2 xlO^'^ at a/s = 1.96 TeV. 

Fig. |3] shows the separate contributions from valence-sea and sea-sea interactions to 
the total differential cross sections for W, Z boson production. For the pp process, the 
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FIG. 3: The rapidity distribution of the differential cross sections for the production of weak 
bosons, W and Z. The upper panels are for ^/s = 7 TeV and the lower panel for ^/s = 14 TeV. 
The contribution from valence-sea (dashed) and sea-sea (dotted) interactions are shown along with 
the total cross section (solid). 

sea-sea interactions for W^ and W^ are the same. As indicated in the figure, the sea-sea 
contributions are dominant in the central rapidity region, while the valence-sea contributions 
dominate in the forward and backward regions. The sea-sea contributions increase with 
collision energy. At rapidity zero the sea-sea process yields 75 % and 86 % of the W^ and 
W^ production cross sections for ^/s = 7 TeV, and 86 % and 92 % for y/s = 14 TeV. 
And we note that because the CKM elements squared have values of 0.95 for |Ktdp and 
|V"csP, and 0.05 for |\4sp and |V"cdP, the valence-sea contributions are dominated by the u^d 
channel for W^ and ud^ for W~ production. For Z production, the sea-sea interactions 
yield 80 % and 89 % for ^/s = 7 and 14 TeV, respectively. Finally we note that the b quark 
interactions have an effect on the Z distributions as large as 4.1 % and 5.8 % in the central 
region, and contribute less than 1 % at \Y\ > 2.5 and |V| > 3.3 for y/s = 7 TeV and 14 TeV, 
respectively. On the other hand, for W production, the b quark contributions are much 
more suppressed because of the small CKM matrix elements squared, with values of 10^^ 
for iKtbp and 2 x 10~^ for |K;6p- Therefore, we can totally ignore the b quark contributions 
for the charge asymmetry of W^ production, and safely ignore for the ratio of W/Z. 
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III. THE CHARGE ASYMMETRY OF W^ PRODUCTION 



In this section, we investigate the rapidity dependence of the charge asymmetry. The 
W^ is produced mainly by the u and d, and W~ by u and d channels. Since the PDFs 
of these quarks are different, a charge asymmetry appears in the W^ and W~ production 
distributions, which is defined as 



Aw{Y) 



daw+ /dY — daw- /dY 



(6) 



da\iY+ /dY + da^- /dY 

The charge asymmetry in pp collision can be approximated in terms of the u and d quark 
distributions as 



AiY) « 
where the weight ui and uj2 are 



Ul 



u{xi) — d{xi) u{x2) — d{x2) \ 



u{xi) + d{xi 



+ UJ2 



U{X2) +d{x2) j 



(7) 



Ui 



002 = 



u(x2)+d(x2) 
u{x2)+d{x2) 



u{xi)+d{xi) . u{x2)+d{x2) 
u{xi)+d{xi) ' u{x2)+d{x2) 

u{xi)+d(xi) 
u(xi)+d{xi) 



u{xi)+d{xi) I u{x2)+d{x2) 
u{xi)+d{xi) ' u{x2)+d{x2) 



(8) 



In Eq. [7] we made the approximation that u{x) = d{x). As noted before, the h quark 
interactions contribute a negligible amount as a result of the very small values of the CKM 
matrix elements^ . The weights cui and UJ2 are shown in Fig. H] as a function of W rapidity. 
Note that ui and uj2 both approach the value of 0.5 when F = 0, i.e. xi = X2- Also, ui 
approaches 1 (0) for large positive (negative) rapidity, while a;2 approaches 1 (0) for large 
negative (positive) rapidity. 



The expression corresponding to Eq. [7] for pp collisions is given by 20] 

u{xi) — d{xi) u{x2) — d{x2) 



AiY) 



D- 



(9) 



_u{xi) + d{xi) u{x2) + d{x2)_ 

The multiplication factor has a value of 0.91 in the central region 20|, in contrast with a 
value of 0.5 for the u factors in pp collisions. From Eq. [9l it follows that A{Y) = —A{—Y) 



^ We can perfectly neglect the c quark contribution to the numerator of A(Y); however, in the denominator 
there are terms from cs and cs which cannot be ignored; they change the normalization. We use the 
original definitions, Eqs. (6) and (14), for our numerical presentation. 
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FIG. 4: The weights ui (dark black) and UJ2 (bright red) as a function of rapidity for ^/s = 7 and 
14 TeV. 

while A(Y) = A{—Y) for pp collisions. In summary, in pp collisions the W~^ distribution 
is symmetric with the W~ distribution, daw±{Y) = daw^{—y)- Therefore, the rate of the 
two is the same aX Y = 0. In contrast, in pp collisions, the rapidity distributions of W~^ 
and W~ are each symmetric with respect to rapidity, daw±{y) = daw±{—y), and their 
magnitudes do not coincide at zero rapidity. Therefore, in the central region for pp collision 

u{x) — d{x) Uy{x) — dj,{x) 



A{Y = 0) 



(10) 



u{x) + d{x) Uy{x) + dv{x) + 2S{x) 
with S{x) = u{x) = d{x), where x = Mw/y/s, while A{Y = 0)=0 in pp collisions. This 
difference makes the determination of u/d possible using pp data. 

Inclusively produced weak bosons decay to leptons, W — )■ lui and Z ^ U, and the ra- 
pidity Yw,z are reconstructed from the energy (Ei) and the longitudinal momentum (pi^l) 
of the secondary leptons. The rapidity of the Z boson can be fully reconstructed from the 
measurement of its decay products. In contrast, the W rapidity cannot be reconstructed be- 
cause the longitudinal momentum of the neutrino is not measured. We therefore reconstruct 
the W rapidity from the momentum of the charged lepton and the (indirectly determined) 
transverse momentum of neutrino [.26]. We introduce the definitions 

1, Ei+pi,L + Pu,,Te^''^ 



Y+ = -In- 



Y,^ ' 



2 El - pi^L + PuuTe ^''± 



(11) 



where 



y;^ = Yi± ln[l + 5 + V'5(2 + 5)], 



(12) 



with 



6 



M^-M^ 



'^Pl,TPui 



T 



(13) 



Here M| is 
{Pl,T + Pu,t] 



the transverse mass of the lepton and neutrino defined as M^ = {\pi^t\ + \Pu,t\Y ■ 
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FIG. 5: The distributions of the reconstructed rapidities Y± and the charged lepton rapidity Yi as 
function of Y. 

Fig. [5] shows the distribution of the reconstructed rapidity Y± and the charged lepton 
rapidity Yj as a function of the rapidity of W. Yi and Y± are related: Yi = (F+ + Y_)/2. 
In the massless limit, the lepton rapidity Yi is the same as the pseudorapidity rj, which is 



23 



25 



34 



42|. As indicated in Fig. O Yi 



used to investigate the lepton charge asymmetry 
correlates roughly with the rapidity Y, therefore, Yi could also be an approximate optional 
choice for Y. In this paper, however, we use the experimentally fully reconstructible Y± to 
investigate the charge asymmetry and the ratio of W/Z. 

The experimentally reconstructed rapidity Y± has a twofold ambiguity because of the 
two solutions for Y^^ from Eq. [121 This can be displayed by comparing the Y± and Y 
distributions: Fig. |6] shows the differential cross sections for W~^ and W~ as a function 
of Y and Y^ at ^/s = 7 TeV and 14 TeV. Thus, we can chose either y+ or Y^ as the 
appropriate observable to study the charge asymmetry. In this paper, we take YI and will 
therefore present results for the positive rapidity range only. In Fig. |6]the solid lines show 
the Y distributions, while the dashed and dotted lines show the Y_ distributions with the 
transverse momentum limited. 
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FIG. 6: The differential cross section for W~^ and W^ as function of the rapidity. The sohd hnes 
indicate the distributions for the real rapidity of the W boson, while the dashed lines are for the 
reconstructed rapidity Y^ with or without the Pt cut. 

We can now redefine the W charge asymmetry A{Y) in terms of the reconstructed rapidity 

Y± as 

daw+/dY_(^^) — daiY-/dY_^^) 



AiY) 



(14) 



daw+dY_(^+) + (ia^y-/rfF_(+) 

for Y > {Y < 0). The resuhs for the asymmetry following Eq. [14] are presented in Fig. 
[7] for different sets of PDFs. For this numerical evaluation, we have included all flavors in 
Eq. El We note that the results for the NLO PDFs (red solid) |29J and NNLO PDFs (blue 
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FIG. 7: The charge asymmetry A{Y^ for the different PDFs for energies of \fs = 7 and 14 TeV as 
a function of the reconstructed rapidity y_. The results for the PDFs are differentiated by color: 



red for MSTW2008-NLO [29], blue for MSTW2008-NNLO 

for CTIO [33]. 



green for CTEQ6.6 



321 ] and black 
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dashed) [29|] are essentially degenerate, as expected. 



IV. CHARM QUARK DISTRIBUTIONS OF THE PROTON 



As shown in sec. II, sea-sea quark contributions significantly contribute to LHC rapidity 
distributions of the weak bosons. We will focus next on the charm contributions cs and cd 
(cs and cd) for W^{W~) and of cc for Z. 





FIG. 8: The charm quark contribution to the rapidity distribution of da/dY for the production of 
W'^,W-, and Z boson with y/s = 1 TeV and 14 TeV. 

The shaded area in Fig. [H] displays the charm quark contributions to the weak boson 
rapidity distributions. As already mentioned, for W production the terms proportional to 
\Vcd\^ are quite small and therefore the dominant subprocess is cs [sc] for W^ iW~)- This 
charm contribution dominates in the central region and contributes 29% to W^ , 33% to W~ , 
and 8% to Z production cross sections at ^/s = 7 TeV. At 14 TeV its infiuence increases to 
36%, 38%, and 11%, respectively. 

In this section, we will investigate the charm PDF using the ratio of the W and Z 
differential rapidity distributions. We redefine the quantity B(Y), introduced in Ref. 

for pp collisions, in terms of the experimentally reconstructed rapidity as 

daw+/dY_ + da\Y-/dY^ 



26| 



BiY) 



daz/dY 



(15) 
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B{Y) is shown in Fig. |9]for different PDFs. As is the case for the charge asymmetry, the 
values of B{Y) evaluated with MSTW2008-NLO PDFs ^ and MSTW2008-NNLO [29\ are 
almost identical. The effect of the cut Pt > 35 is shown by the dashed lines. The values of 
B(Y) are 3.05 - 3.11 for the CTEQ6.6 PDFs ^, and 3.07- 3.10 for MSTW2008 for |r| < 2 
at y/s = 7 TeV. At v^i = 14 TeV, they are slightly reduced, to 2.96 - 3.11 for CTEQ PDFs, 
and 2.99 - 3.09 for MSTW PDFs in the range |y| < 3. The results with the transverse 
momentum Pt applied differ by 1.3% - 4.3% for |r| < 2 at v^i = 7 TeV, and 1.3% - 4.4% 
for |F| < 3 at v^ = 14 TeV. 

In Fig. [TOl B(Y), given by Eq. [151 is shown as a function of the charm contribution e, 
which is defined as 



c(x) 



2c(x) 



{{u{x) + d{x))/2) {u{x)+d{x)) 



(16) 



No cut on the transverse momentum pt is applied. Because the NNLO set gives almost 
the same results as the NLO for MSTW PDFs [29|, we only show the results for the NLO 
PDFs. Comparing the results with and without the contributions of charm we find that e 
has a value of 0.5 -0.75. 

In addition to B{Y), we calculate the the ratio of the integrated cross sections, B{e), 
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FIG. 9: B{Y^ as a function of rapidity for ^fs = 7 and 14 TeV. The distributions without Pt cut 
are displayed by solid lines while the results evaluated with Pt > 35 GeV are shown asdashed 



lines. As in Fig. [3 the CTEQ6.6 PDF 
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FIG. 10: BiJ) as a function of e (e = 2c(x)/(n(x) + d{x))) for ^ = 7 and 14 TeV. The CTEQ6.6 
and MSTW 2008-NLO PDF sets are used. The red sohd hues show the results including the charm 
structure function in the sample PDFs. 



which is given by 



5(e) 



doy^Jr jdY^ + doyy- jdY^ 



(17) 



vith y cut doz/dY 

The result is displayed in Fig. [TT]as function of e. B{t) is evaluated for |y| = 1 and 3 for 
^/s = 7 TeV, and 2 and 4 for 14 TeV. The crosses in the figure indicate the B{e) values 
evaluated with the charm distributions from the PDF sets. For example, B{e) evaluated for 
|r| < 3 is 3.07 for CTEQ6.6 [3^ and 3.08 for MSTW2008 PDFs y at ^i = 7 TeV. At 
y/s = U TeV, we find 3.06 for both PDFs for |r| < 4. These values yield e 0.63 to 0.60 
at yi = 7 TeV, and 0.68 to 0.66 at 14 TeV for the CTEQ6.6 and MSTW2008 PDF sets, 
respectively. 
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FIG. 11: The ratio of the integrated cross sections B{e) as a function of e for different ranges of Y. 
The values of B{e) evaluated with the charm distributions from each sample PDF are represented 
by the crosses. 

V. SUMMARY AND CONCLUSIONS 



Using sample PDFs, we have examined the charge asymmetry for u and d quarks as 
well as the ratio of cross sections B{Y) and Bie) in the experimentally accessible weak 
boson rapidity range at the LHC We subsequently focused on the role of charm quark 
contributions. For the charge asymmetry the effects of the heavy quark PDFs c(x) and h{x) 
are negligible, even at the highest energy. One can therefore isolate and explore the u[x) 
and d{x) (or u.u{x) and du{x)) structure functions. In this context, we presented a simphfied 
expression of the charge asymmetry in terms of the u and d distribution functions. 

We also studied the charm contribution e = 2c(x)/('u(x) + d{£)) to B{Y) and B(t). The 
CTEQ6.6 M and MSTW2008 PDF [29| sets indicate a value of e of about 0.6 for ^s = 
7 TeV, and ~ 0.7 for 14 TeV. The measurement of the rapidity distributions of W/Z (or 
B{Y)) interpreted in terms of B{Y) and B{e) will lead to a straightforward determination 
of the charm quark component of the the sea quarks at the LHC. 
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